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I. INTRODUCTION 

Microscopic theories for calculating nuclear masses 
and/or binding energies (see, e.g., [1-3]), have been re- 
vived and further elaborated with the advance of com- 
putational resources. These advances are now sufficient 
to perform global studies based on, e.g., self-consistent 
mean field theory, sometimes also denoted by density 
functional theory (DFT) [4, 5]. One particular aspect 
of the nuclear binding problem is a phenomenon of odd- 
even staggering (OES) of the binding energy. Numerous 
microscopic calculations have been published that treat 
individual isotope chains. However, it might be necessary 
to examine the whole body of OES data to draw general 
conclusions [6]. 

Theoretically, OES values are often inferred from the 
average HFBCS or Hartree-Fock-Bogoliubov gaps [7-9], 
rather than directly calculated from the experimental 
binding energy differences between even and odd nuclei. 
It should be mentioned that the average HFB gaps are 
sometimes substantially different from the odd-even mass 
differences calculated from experimental binding ener- 
gies. In this work, we compare directly the calculated 
OES with the ones extracted from experiment. One 
should say that there are also several prescriptions to ob- 
tain the OES from experiments, such as 3-point, 4-point, 
and 5-point formulas [6]. We adopt the 3-point formula 
A*^^' centered at an odd nucleus, i.e., odd-N nucleus for 
neutron gap and odd-Z nucleus for proton gap [2]: 



A(3)(Ar, z) 



B{N - l,Z) - 2B{N,Z) (1) 
+ B(7V + 1,Z)1 , 



where B{N, Z) is the binding energy of {N, Z) nucleus 
and IT A = {^)^ is the number parity with A = N + Z . 

For even nuclei, the OES is known to be sensitive not 
only to the pairing gap, but also to mean field effects, 
i.e., shell effects and deformations [6, 7]. Therefore, the 
comparison of a theoretical pairing gap with OES should 



be done with some discretion. One advantage of AJ, 
(iV = odd in Eq. (1)) is the suppression of the contri- 
butions from the mean field to the gap energy. Another 
advantage of /sSf^ [N, Z) is that it can be applied to more 
experimental mass data than the higher order OES for- 
mulas. At a shell closure, the OES (Eq. (1)) does not go 
to zero as expected, but it increases substantially. This 
large gap is an artifact due to the shell effect, which is 
totally independent of the pairing gap itself. 

Recently, an effective isospin dependent pairing inter- 
action was proposed from the study of nuclear matter 
pairing gaps calculated by realistic nucleon-nucleon in- 
teractions. In Ref. [8], the density— dependent pairing 
interaction was defined as 



2) = Vo gr[p, /3t,] ,5(ri - v^), 



(2) 



where p = pn -\- pp \s the nuclear density and /3 is the 
asymmetry parameter /3 — {pn — Pp)/ P- The isovector 
dependence is introduced through the density-dependent 
term g^. The function g,- is determined by the pairing 
gaps in nuclear matter and its functional form is given 

by 



lr[p.PT,] = l-UpT,)T^, 



,(3) 



where pQ—0.16 fm is the saturation density of symmet- 
ric nuclear matter. We choose isi^Tz) = 1 — ini^Tz) and 
fn(/3r^) = = [pn(r) - pp{r)] r^/pir). The parameters 
for gr are obtained from the fit to the pairing gaps in 
symmetric and neutron matter obtained by the micro- 
scopic nucleon-nucleon interaction. 

In the literature and in many mean field codes publicly 
available such as the original EV8 code [10], a pure con- 
tact interaction is used without an isospin dependence. 
In our notation, this amounts replacing the isospin de- 
pendent function g^ in Eq. (2) by the isoscalar function 



= 1 



'7s 1 — 



(4) 
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TABLE I: Parameters for the density-dependent function gr defined in Eqs. (2) and (3) for the IS+IV interaction (first row) 
and gs in Eq. (4) for the IS interaction. The parameters for gi- are obtained from the fit to the pairing gaps in symmetric and 
neutron matter obtained with the microscopic nucleon-nucleon interaction. The paring strength Vb is adjusted to give the best 
fit to odd-even staggering of nuclear masses. The parameters for gs correspond to a surface peaked pairing interaction with no 
isospin dependence. 



interaction Vq (MeVfm^) po (fm ^) Vb cts Vn cin 

gr (isotopes) 1040 0.16 0.677 0.365 0.931 0.378 

g^ (isotones) 1120 0.16 0.677 0.365 0.931 0.378 

gs (isotopes) 1300 0.16 1. 1. — — 

gs (isotones) 1500 0.16 1. 1. — — 
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FIG. 1: (Color online) The mean square deviation a of OES 
between experimental data and the HF-l-BCS calculations. 
The filled circles and squares correspond to the results with 
IS pairing for neutron and proton gaps, respectively, while the 
filled diamonds and triangles are those of IS-I-IV pairing for 
neutron and proton gaps. Experimental data are taken from 
Ref. [17]. See the text for details. 



The EV8 code has been modified, using the filling ap- 
proximation, to account for mass calculations for odd-N 
and odd-Z nuclei and is publicly available as the EVSodd 
code [11]. It has also been modified to include isospin de- 
pendent paring, by means of Eq. (2). The parameters 
of the isoscalar interaction were adjusted with EV8 to a 
best global fit of nuclear masses [12]. They correspond 
to a surface peaked pairing interaction (Eq. (4) with r]s 
not too far from the unity). 

A recent publication has explored the isospin depen- 
dence of the pairing force for the OES effect for a few 
selected isotopic and isotonic chains [13]. Here we have 
made a more ambitious study by extending the calcula- 
tion to the whole nuclear chart. We have also explored 



several other observables such as neutron and proton 
radii systematically which may allow for more solid con- 
clusions on isospin dependent pairing interactions. 

This paper is organized as follows. In section II we dis- 
cuss our numerical calculation strategy. Our results are 
presented in section III for the energies, separation en- 
ergies, OES energies, and nuclear radii. Our conclusions 
are presented in section IV. 

II. CALCULATION STRATEGY 

The HF-l-BCS calculations are performed by using 
SLy4 Skyrme interaction which was found to be the most 
accurate interaction for studying OES for a few selected 
{N = 50, 82) isotonic and (Sn and Pb)isotopic chains [13]. 
Our iteration procedure used in connection to EV8odd 
achieves an accuracy of about 100 keV, or less, with 500 
Hartree-Fock iterations for each nuclear state. Our cal- 
culations were performed with the now decommissioned 
XT4 Jaguar supercomputer at ORNL, as part of the 
UNEDF-SciDAC-2 collaboration [14]. 

The HF-l-BCS calculations were first performed for 
even-even nuclei. The variables in the theory are the or- 
bital wa ve func tions (f>i and the BCS amplitudes Vi and 
Ui — y^l — vf. By solving the BCS equations for the 
amplitudes, one obtains the pairing energy from 

Epair = ^ VijUiViU-jVj + ^ Vuvf (5) 

where Vij are the matrix elements of the pairing interac- 
tion, Eq. (2), namely 

V.j^Vo J d'r\Mr)f\cl>j{r)fgr[p{r),P{r)T,], 

where p(r) = ^'^^l</'i(^)P■ 
After determining the single-particle energies of even- 
even nuclei, the odd-A nuclei are calculated with the so- 
called filling approximation for the odd particle starting 
from the HF-l-BCS solutions of neighboring even-even nu- 
clei: ones selects a pair of i and i orbitals to be blocked. 
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FIG. 2: (Color online) Binding energy differences between experimental data and calculations using the HF+BCS model with 
the IS and IS+IV pairing interactions. The left panels show the differences for even Z isotopes varying neutron numbers 
including both odd and even numbers. The right panel show those for even N isotones varying proton numbers including both 
odd and even numbers. The thin lines show the closed shells at N(Z) = 20, 28, 40, 50, 64, 82 and 126. Experimental data are 
taken from Ref. [17]. See the text for details. 



and chanecs the BCS parameters vf and for these 

'' i 

orbitals. The change is to set ~ ^ ^/"^ 
(5) for the pairing energy at an orbital near the Fermi 
energy. Note that this approximation gives equal occu- 
pation numbers to both time-reversed partners, and does 
not account for the effects of time-odd fields. More de- 
tails of the procedure are presented in Ref. [12]. 

The effect of time-odd HF fields on the mass were stud- 
ied in Refs. [15, 16]. It was pointed out that the effect 
of the time-odd fields is of the order of 100 keV for the 
binding energy depending strongly on the configuration 
of the last particle, and does not show any clear sign of 
isospin dependence. Thus the time-odd field might not 
change conclusions of the present study in the following, 
while quantitative accuracy might need some fine tuning 
of the pairing parameters. 

For the pairing channels we have taken the surface- 
type contact interaction, Eq. (4), and the isospin de- 
pendent interaction, Eq. (3). The density dependence 
of the latter one is essentially the mixed-type interaction 
between the surface and the volume types. The pairing 
strength Vq depends on the energy window adopted for 
BCS calculations. The odd nucleus is treated in the fill- 
ing approximation, by blocking one of the orbitals. The 
blocking candidates are chosen within an energy window 
of 10 MeV around the Fermi energy. This energy win- 
dow is rather small, but it is the maximum allowed by 
the program EVSodd. It is shown that the BCS model 
used in the EVSodd code gives almost equivalent results 



to the HF-I-Bogoliubov model with a larger energy win- 
dow, except for unstable nuclei very close to the neutron 
drip line [12]. The pairing strengths Vq for IS and IS-flV 
pairing interactions are adjusted to give the best fit to 
odd-even staggering of nuclear masses in a wide region of 
the mass table. 



III. NUMERICAL RESULTS 
A. Global data on odd-even staggering 

The results for the mean square deviation of our global 
mass table calculations are shown in Fig. 1. Table I gives 
the values Vq in Eq. (2) and the parameters for g^ and 
gs is Eqs. (3) and (4) used in the present work. Optimal 
pairing strength values were found to be different for iso- 
tones (varying Z , constant N) and for isotopes (varying 
iV, constant Z). 

Figure 1 shows the mean square deviation a of OES be- 
tween experimental data and the HF-I-BCS calculations. 
The mean square deviation a is defined as 



\ 



Ni 



Af ^ [HF + BCS) - Af ^ [exp] /N, (6) 



where Ni is the number of data points. For the IS interac- 
tion, the results for neutrons show a shallow minimum at 
Vq ^ (1100 — 1300) MeV-fm'^. For protons, the minimum 
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FIG. 3: (Color online) The same as Fig. 2 but for neutron and proton separation energies. See the caption to Fig. 2 and the 
text for details. 



becomes at around Vq ^ 1500 MeV-fm''. This clifFerence 
makes it difficult to determine a unique pairing strength 
common for both neutrons and protons. The results of 
IS+IV pairing show a minimum at Vq ~ 1100 MeV-fm'^ 
for both neutron and proton OES which makes it easier 
to determine the value for the pairing strength. Adopted 
values for the following calculations are listed in Table I. 

The systematic study of HF+BCS calculations are per- 
formed for various isotopes and isotoncs for all available 
data sets with = 8, • • • , 102) and (A^ = 8, • • • , 156), 
respectively. Binding energy differences between experi- 
mental data and HF+BCS 

SE = lE'^^'P ~ £""^1 (7) 

are shown for both IS and IS+IV pairing interactions 
in Fig. 2. The left panels show the values SE varying 
neutron numbers (including both odd and even numbers) 
for each even Z. The right panels show the values SE 
varying proton numbers (including both odd and even 
numbers) for each even N. With IS pairing, we can see a 
rather large deviation for Z = 50 isotopes in the upper 
right panel. This difference disappears in the case of 
IS+IV pairing shown in the lower right panel. On the 
other hand, for the A = 82 nuclei, the IS+IV interaction 
does not work that well. As far as the binding energies 
are concerned, the best results with IS+IV interaction are 
obtained for nuclei with iV = 60 - 78 and A^ = 86 - 96. 

Separation energy differences between experimental 
data and HF+BCS for protons and neutrons are plotted 
in Fig. 3. The HF+BCS results of neutron separation 
energies Sn are reasonable for medium and heavy mass 
nuclei with A^ = 60 — 120, except near the closed shell 



A = 82. For heavy nuclei with A = 126, the calculated 
results are poorer than in other mass regions. For proton 
separation energy Sp, the HF+BCS also gives reasonable 
results, except in the Z — 50 and 82 mass regions. 

In order to see the different outcomes between IS and 
IS+IV pairing interactions, the HF+BCS model calcula- 
tions are shown together with empirical data in Fig. 4. 
In most of cases, the difference between the two pairing 
interactions are small. However, we can see a clear im- 
provement of the agreement of Sp with empirical data of 
A^ = 136 isotones with IS+IV pairing in Fig. 5. 

The differences of neutron OES A„ and proton OES 
Ap between HF+BCS and empirical data are shown in 
Fig. 6 for both IS and IS+IV pairing, respectively, in 
the upper and lower panels. The agreement between 
HF+BCS and the empirical data are good in the overall 
mass region except for masses with Z = 50 and at a small 
mass region A < 60. To clarify the difference between IS 
and IS+IV pairing, the OES differences A„ are shown in 
the upper panel of Fig. 7 for Z = 52, 78 and 92 isotopes. 
The HF+BCS results are compared with the experimen- 
tal data and also the phenomenological parameterization 
based on liquid drop model, 

A = c/A" (8) 

with c = 4.66(4.31) MeV for neutrons (protons) and 
a = 0.31 which gives the rms residual of 0.25 MeV [12]. 
We can see clearly a better agreement of IS+IV results 
with empirical data for all isotopes. In the lower panel 
of Fig. 7, the OES differences A„ are shown for A = 76, 
102 and 112 isotones. The results with IS+IV pairing 
certainly improve systematically the agreement with em- 
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FIG. 4: (Color online) Neutron separation 6I16rgi6S Sn of three isotope chains with Z = 46, 72 and 90 calculated by IS and 
IS+IV interactions in HF+BCS model (upper panels). The lower panels show the proton separation energies Sp for A'' = 50, 
90 and 136 isotones. Experimental data are taken from Ref. [17]. See the text for details. 



pirical data, especially for N — 102 isotones. The large 
increase of the HF+BCS model results at Z = 81 is an 
artifact due to the shell closure a.t Z = 82. It is interest- 
ing to notice that the liquid drop formula gives smooth 
mass number dependence which reflects well that of very 
heavy isotones with TV = 112. 

The average gaps are tabulated for high and low 
isospins in Table II. Each isotope (isotone) in Fig. 7 
is divided into two subsets of almost equal numbers of 
nuclei by a cut at some value of / = {N — Z)/A. Both the 
average proton and neutron A*^'^-' show smaller values for 
higher isospin so that the pairing interaction is weaker for 
neutron-rich nuclei. The IS-I-IV interaction reproduces 
properly the difference of the neutron A'^'^' between high 
and low isospin nuclei. For proton A^'^^ also, the IS-I-IV 
pairing gives a good account of the isospin effect than the 
IS pairing. 



B. Nuclear radii 

Nuclear radii provide basic and important information 
for various aspect of nuclear structure problems. The 
proton radii, or equivalently the charge radii with the 
correction of finite proton size, can be determined accu- 
rately by electron scattering and muon scattering exper- 
iments. However it is difficult to determine the neutron 



radii of finite nuclei with the same accuracy level as that 
of the proton radii while there were several experimental 
attempts to determine the difference of the neutron to 
proton radius [19-21]. It should be noticed that the dif- 
ference of the neutron and proton radii, (5r„p = r„ — Tp, is 
called the neutron skin. It is thought that <5r„p can give 
important constrains on the effective interactions used in 
nuclear structure study [22]. 

The neutron and proton radii of various isotopes and 
isotones are calculated by using the HF-I-BCS model 
with the two pairing interactions, IS and IS-I-IV. The 
results of neutron radii are shown in the left panel of 
Fig. 8. Since we do not find any appreciable differences 
between the two pairing interactions in the results, the 
results of IS-I-IV interactions will be mainly discussed 
hereafter. The results obtained by a simple empirical 
formula r„ = roiV^/'^ with tq = 1.139 fm [18] are also 
plotted in the figure. In general, the simple formula for 
r„ agrees well with the HF-f-BCS results. It is noticed 
that the HF-I-BCS model gives larger neutron radii for 
nuclei with < 40 than the simple formula but smaller 
for nuclei with TV > 120. The proton radii for N = 20, 
28, 40, 50, 82 and 126 isotones are shown as a func- 
tion of proton number Z in the right panel. The sim- 
ple Z^/^ dependence is also plotted to follow the formula 
Tp = 1.263/Z^'''^. The simple formula in general gives a 
good account of the HF-I-BCS data and could be a good 
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TABLE II: Average A^^^' for low isospin and high isospin nuclei and its difference. See the text for details. 





Data set 




Low isospin 


High isospin 


Difference 


Neutrons 


Z = 52 


Exp 


1.36 


1.08 


-0.28 






IS 


1.52 


1.41 


-0.11 






IS+IV 


1.40 


1.19 


-0.21 




Z = 78 


Exp 


1.13 


0.99 


-0.14 






IS 


0.96 


1.16 


0.20 






IS+IV 


0.87 


0.91 


0.04 




Z = 92 


Exp 


0.77 


0.56 


-0.21 






IS 


0.90 


0.80 


-0.10 






IS+IV 


0.70 


0.55 


-0.15 


Protons 


N = 76 


Exp 


1.19 


0.93 


-0.26 






IS 


1.13 


0.87 


-0.26 






IS+iV 


1.13 


0.98 


-0.15 




N = 102 


Exp 


0.96 


0.63 


-0.33 






IS 


0.79 


0.39 


-0.40 






IS+IV 


0.92 


0.59 


-0.33 




Z = 112 


Exp 


0.87 


0.66 


-0.21 






IS 


0.58 


0.61 


0.03 






IS+IV 


0.67 


0.70 


0.03 
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FIG. 5: (Color online) Proton separation energies of 
N = 136 isotones calculated by IS and IS+IV interactions 
in HF+BCS model. Experimental data are taken from Ref. 
[17]. See the text for details. 



in heavier isotopes, larger than Z = 56. The results of 
IS and IS+IV pairings are shown in the left panel and 
right panel, respectively. In general, the two pairing in- 
teractions give almost the same results as shown in Fig. 
9. However, it is noticed that the IS+IV pairing gives 
somewhat smaller neutron skins than the IS pairing in 
very neutron-rich nuclei such as ^'^^Sn, ^^°Ba and ^^®Po. 
The calculated values are compared with empirical data 
of Sn isotopes obtained from studies of spin-dipole reso- 
nances [19] and antiprotonic atoms [20] in the left panel 
Fig. 10. The calculated values show reasonable agree- 
ment with the empirical data within the experimental 
error. The isospin dependence of neutron skin is shown 
in the right panel of Fig. 10 together with empirical 
values obtained by antiprotonic atom experiments in a 
wide range of nuclei from '"^Ca to ^'^^U. The slope of ex- 
perimental data as a function of the isospin parameter 
/ = (TV — Z)/A is reproduced well by our calculations. 



starting point for describing the isospin dependence of 
nuclear charge radii. However we can see some deviation 
between the HF+BCS and the simple formula especially 
heavy = 50 and A^ = 82 isotones. 

The neutron skin 7-„ — calculated by HF+BCS model 
with the two pairing interactions are shown in Fig. 9. 
The neutron skin becomes as large as 0.4 fm near the 
neutron drip hne with Z < 28. On the other hand, the 
neutron skin is at most 0.25 fm in neutron-rich nuclei 
with Z > 50. For proton-rich nuclei, the proton skin 
becomes 0.1 fm with Z < 56 and smaller than 0.05 fm 



The neutron skin of ^''^Pb has been discussed inten- 
sively in relation with neutron matter properties. The 
systematic studies of scattering data yield the empirical 
value Tn — Tp = 0.17 ± 0.02 fm which is close to another 
empirical value r„ — Vp = 0.15 ± 0.02 fm from the study 
of antiprotonic-atom systems. The model independent 
determination of parity violation experiment at Jefferson 
Laboratory [22] has been proposed and performed re- 
cently to obtain the neutron skin of ^°®Pb. However the 
statistics was poor and needs improvement by more data 
accumulation. Our calculated value r„ — fp = 0.157 fm 
is close to the experimental values by the two systematic 
studies. 
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FIG. 6: (Color online) The same as in Fig. 2 but for OES for neutrons and protons. See the caption to Fig. 2 and the text for 
details. 



IV. SUMMARY AND CONCLUSIONS 

In summary, we studied the binding energies, separa- 
tion energies and OES by using HF+BCS model with 
SLy4 interactions together with the isospin dependence 
pairing (IS+IV pairing) and isoscalar (IS pairing) in- 
teractions. The calculations are performed with the 
EVSodd code for even-even nuclei and also even-odd 
nuclei using the filling approximation. For the neu- 
tron pairing gaps, the IS-I-IV pairing strength decreases 
gradually as a function of the asymmetry parameter 
(p„(r) — pp(r))/p(r). On the other hand, the pairing 
strength for protons increases for larger values of the 
asymmetry parameter because of the isospin factor in 
Eq. (3). The empirical isotope dependence of the neu- 

tron OES, A„ , is well reproduced by the present calcula- 
tions with the isospin dependent pairing compared with 
the IS pairing. We can also obtain a good agreement 
between the experimental proton OES and the calcula- 
tions with the isospin dependent pairing for = 50 and 
= 82 isotones. 

The neutron and proton radii were also studied by us- 
ing the same HF-I-BCS model with the two pairing inter- 
actions. The two pairing interactions give essentially the 
same results for the radii except for a few very neutron- 
rich nuclei. We found systematically large neutron skins 



in very neutron-rich nuclei with |r„ — rp\ ^ 0.4 fm, while 
the proton skin is rather small even in nuclei close to 
the proton drip line because of the Coulomb interaction. 
The calculated results of neutron skin show reasonable 
agreement with the empirical data including the {N — Z) 
dependence of the data. 

We tested the IS-I-IV pairing for the Skyrme interac- 
tion SLy4 and found the results reproduce well the sys- 
tematical experimental data. Thus, we confirm a clear 
manifestation of the isospin dependence of the pairing 
interaction in the OES in comparison with the experi- 
mental data both for protons and neutrons. 
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FIG. 7: (Color online) The neutron and proton OES, A„ and Ap, calculated with the HF+BCS model with IS and IS+IV 
interactions. Experimental data are taken from Ref. [17]. See the text for details. 
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FIG. 8: (Color online) Neutron and proton radii of various isotopes and isotones calculated by means of the IfF+BCS model 
with IS+IV interaction. The solid lines are empirical fits used in Ref. [18]. See the text for details. 




FIG. 9: (Color online) Neutron skin neutron for isotopes and isotones calculated by means of the HF+BCS model with IS and 
IS+IV interactions. The thin lines indicate the closed shells with N (or Z)=20, 28, 40, 50 ,64, 82 and 126. See the text for 
details. 
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FIG. 10: (Color online) Left - Neutron skin for Sn isotopes obtained with the HF+BCS model with IS+IV interaction. 
Experimental data are taken from Ref. [19, 20]. Right - Neutron skin as a function of isospin parameter / — {N — Z)/A 
calculated by means of the HF+BCS model with IS+IV interaction. Experimental data are taken from Ref. [21]. See the text 
for details. 



